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Abstract 


In  this  paper  we  consider  the  problem  of  locating  a  single 
facility  on  a  tree  graph  so  as  to  minimize  the  maximum  "round  trip" 
distance  from  this  facility  to  a  given  set  of  node  pairs.  The  node 
pairs  might  represent  a  neighborhood  and  a  hospital  or  a  neighborhood 
and  a  police  station  to  which  service  must  be  provided.  Each  node  pair 
can  be  assigned  a  weight  which  represents  the  relative  amount  of  service 
which  must  be  provided  at  chat  node  pair  to  Che  total  amount  of  estimated 
service.  In  Che  case  where  all  of  the  weights  are  equal  we  provide  the 
results  necessary  to  justify  a  simple  0(m)  procedure  for  finding  an 
optimal  location,  where  m  is  the  number  of  node  pairs.  When  the  weights 
are  unequal  an  0((tt+mHgn  +  migm)  procedure  for  finding  an  optimal  location 
is  presented  where  n  is  the  number  of  nodes  in  the  tree  graph.  The 
classical  one  center  problem  on  a  tree  graph  is  a  special  case  of  this 
problem.  Examples  and  references  to  related  research  are  included.. 


1.  Problem  Formulaelon 

Given  a  connected  undirected  network,  N,  with  node  or  vertex  set 
V,  and  edge  or  arc  set  E,  let  d(x,y)  be  the  length  of  a  shortest 
path  between  any  two  points  x  and  y  on  N  (x  and  y  are  not  necessarily 


nodes ) 

.  Let  {(p.^ 

> 1 P±  €  V,  qt  6  V,  1  < 

i  <  m} 

be  a  given  set 

of  node 

pairs . 

For  any  point 

x  €  N  we  define. 

rt(x)  -  d(x,pt)  +  d(x,qi)  +  ^ 

1  <  i 

<  m 

(1) 

where 

r . (x)  is  the 

i 

"round  trip"  distance 

from  x 

to  node  pair 

i,  and 

k^  is 

a  known  constant.  In  particular  if 

we  set 

kt  -  d(pi,qi) 

then 

rt(x) 

represents  the 

length  of  a  shortest 

path  from  x  to  p^, 

from 

to  q^,  and  then  from  q^  back  to  x,  thus  the  name  round  trip.  In 
the  context  of  this  paper  we  let  k^  -  ^(p^.q^)  +  k^,  so  that  r^(x) 
represents  the  round  trip  distance  from  x  to  node  pair  i  plus  a 
fixed  charge. 

Now  suppose  we  associate  with  each  node  pair,  (p^.q^),  a  positive 
weight  w^  which  denotes  the  importance  or  relative  value  of  a  round 
trip  to  node  pair  i.  Next  we  define, 

r (x)  -  max  [w.r  (x)]  (2) 

1  <  i  <  a  11 

where  r(x)  is  the  maximum  weighted  round  trip  distance  from  x  to 
any  node  pair.  Now  we  consider  the  problem, 

P„:  Minimize  r(x)  .  (3) 

N  x  €  N 

In  problem  we  seek  the  set  of  locations,  RC,  on  N  such  that 

the  maximum  weighted  round  trip  distance  from  any  point  in  RC  is 
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minimized  over  all  possible  locations  on  N.  In  this  paper  we  address 
the  special  case  of  problem  where  N  is  a  tree,  that  is  we  consider, 

P  :  Minimize  r(x)  (4) 

T  x  6  T 

where  T  is  a  connected  graph  with  no  cycles. 

Applications  of  round  trip  location  can  be  found  in  many  situations, 
such  as  in  public  sector  location  problems,  where  locating  a  facility 
with  the  objective  of  minimizing  the  worst  case  occurrence  is  a  valid 
criterion.  For  example  consider  the  location  of  an  ambulance  station  in 
an  urban  center  where  the  edges  or  arcs  of  the  network  correspond  to 
streets  and  highways,  the  nodes,  p^,  represent  neighborhoods  where 
the  demand  for  ambulance  service  will  occur  and  the  nodes,  q^,  cor¬ 
respond  to  local  hospitals  or  service  centers.  In  general  a  patient 
picked  up  by  an  ambulance  in  neighborhood  p^  would  be  transported  to 
the  nearest  hospital  and  thus  node  p^  would  be  paired  with  a  node  q.^ 
representing  this  nearest  hospital.  However,  in  some  cases  the  resources 
necessary  to  treat  a  patient,  such  as  equipment,  specialists,  or  drugs, 
are  unavailable  at  the  nearest  hospital  (this  is  especially  true  in  the 
case  where  local  hospitals  share  resources)  and  thus  it  might  be  necessary 
to  transport  the  patient  to  a  hospital  which  is  further  away.  Thus  a 
neighborhood  might  be  paired  with  several  hospitals  besides  the  one 
nearest  to  it.  In  order  to  distinguish  the  relative  frequency  of  occur¬ 
rence  or  importance  of  each  neighborhood-hospital  pair  we  assign  weights 
to  each  pair,  the  larger  weights  being  given  to  the  neighborhood-hospital 
pairs  with  the  highest  priority.  Then  the  objective  in  P^  is  to  locate 
an  ambulance  station  so  that  the  maximum  weighted  round  trip  distance 


-3- 


from  chis  ambulance  station  to  a  neighborhood-hospital  pair  is  minimized. 

For  other  discussions  of  round  trip  location  see  [1]  and  [2]. 

In  [ 1]  ,  Chan  and  Francis  described  an  algorithm  of  complexity 

2 

0(m  )  for  finding  an  optimal  solution  to  P^.  In  section  two  we 
provide  a  brief  summary  of  their  approach.  Then  in  section  three  we 
provide  results  which  show  that  for  an  unweighted  tree,  a  very  simple 
procedure  with  complexity  0(m)  can  be  used  to  find  an  optimal  solution 
to  P^.  In  section  four  we  provide  a  method  for  finding  an  optimal  solu¬ 
tion  to  the  weighted  version  of  PT  in  0((n+m)-t.gn  +  migm)  (where  n  is  the 
number  of  nodes  in  T)  which  is  based  on  the  concept  of  a  "local1'  round 
trip  center  and  the  "centroid"  of  a  tree. 

Finally  we  would  like  to  point  out  that  in  the  special  case  of  P^ 
where  pt  ■  and  k^  =■  0  for  all  i,  i  <  i  <  m,  P^  reduces  to 

the  well  known  one  center  problem  which  has  attracted  a  considerable 
amount  of  attention  in  the  literature.  In  particular  part  of  the  work 
in  [ 6]  and  [ 7]  provided  strong  motivation  for  the  research  presented 
in  sections  three  and  four  of  this  paper  respectively. 

2 .  Notation  and  Summary  of  Previous  Results 

In  this  section  we  introduce  some  of  the  notation  and  definitions 
which  are  used  throughout  the  remaining  sections  of  this  paper  plus  a 
brief  summary  of  some  previous  results  for  P^  which  can  be  found  in 

U1 . 

Notation: 

T  -  a  tree;  a  connected  graph  with  no  cycles. 

V  -  the  set  of  nodes  of  T;  | vj  *  n 

[(Pl^i)}  -  a  set  of  given  node  pairs  1  <  i  <  m 
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(Note: 

Pt(x) 

Ptj 

(Not# : 

P(x,y) 

V(x) 

Tv(x) 


RC 


-  Che  (unique)  path  connecting  p,  and 

-  (z  €  Tl(d(pi>z)  +  d(z,qi)  ■  d (pi , qt ) } 

•  The  point  on  P^  which  is  closest  to  x. 

•  che  point  on  P^  which  is  closest  to  every  point  on  p^ 

If  P ^  H  Fj  ^  ®  then  let  p^  ■  p^  be  the  point  on  the  connected 
path  P^  H  P  which  is  closest  to  p^.) 
the  (unique)  path  connecting  x  and  y 

-  che  sec  of  all  nodes  in  V  which  are  adjacent  to  x 

-  defined  for  all  v  €  V(x),  is  the  connected  subtree  of  T 

which  contains  che  segment  (x,v)  plus  all  the  points  in  T 
which  are  connected  to  x  by  a  path  containing  v. 

•  che  sec  of  all  optimal  solutions  to  P_. 


Sow  we  give  a  brief  outline  of  a  procedure  for  calculating  RC,  the 
sec  of  optimal  solutions  to  P^-  First  some  definitions  are  necessary. 
We  define 

(d(pt.qt)  +  ^  ) 

gt  -  - j -  1  <  i  <  m  (5) 

where  g^  represents  one  half  of  che  fixed  charge  portion  of  che  round 
crip  distance  to  node  pair  i. 

Next  we  note  that  for  all  x  €  T  and  for  all  y  €  p^  that 
p^(x)  €  P(y,x).  In  particular  if  x  €  p^  then  p^(x)  »  x.  Given  this 
property  and  (5)  one  can  easily  verify  that  r^(x)  can  be  redefined 


from  (1)  as 
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ri(x)  -  2(d(x,pi(x>)  +  gj  .  (6) 

This  alternate  definition  of  r^(x)  will  be  used  in  subsequent 
sections . 

Now  we  define  a  number ,  H,  which  is  shown  in  [1]  to  equal 
the  optimal  objective  function  value  for 
For  1  <  i  <  j  <  m  define 


h.  . 
ij 


2(w1w,P,i  +  w.g.  +w.gt) 


(w .  +  w .  ) 
1  J 


(7) 


and 


H  =»  h_  *■  max  £ h  }  (8) 

Note:  for  i  =  j,  h^  ■  2g^. 

Intuitively  H  represents  the  weighted  round  trip  "diameter" 
of  the  tree.  Since  H  equals  the  optimal  objective  function  value 
for  PT>  the  set  of  optimal  solutions,  RC,  can  be  written  as 

RC  *  [x|r. (x)  <  —  ,  1  <  i  <  m]  (9) 

1.  ”  w . 

1 

(H  -  2w.g  ) 

or  letting  u^^  *  - . .  1  and  considering  the  definition  of 

r^(x)  in  (6)  we  have 

RC  *  £x {  d(x,pi(x)>  <  u^  i  <  i  <  m]  •  (10) 

Remark  1:  It  is  well  known  (see  [3] )  that  minimizing  a  convex  function, 

such  as  r(x),  on  a  tree  yields  a  set  of  optimal  solutions  which  is 
convex.  By  convex  on  a  tree  we  mean  for  any 
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€  T  let  z  be  the  (unique)  point  which  is 
\-d(x^,x^)  unit3  from  x^  and  (l-\ ) -d (x^.Xj)  units 
from  *2’  0  <  X  <  1;  then  T  is  convex  iff  z  €  Pfa^.x^). 

Remark  2;  Convexity  and  connectedness  are  equivalent  concepts  when 
defined  on  a  tree. 

From  Remarks  1  and  2  we  see  that  RC  is  a  convex,  or  equivalently 
a  connected  sec  on  T. 

Now  the  method  described  in  [1]  to  find  RC  works  as  follows: 

(i)  first  compute  H,  and  (ii)  then  compute  RC  depending  upon  which 
of  the  following  two  mutually  exclusive  cases  occur, 

Case  1:  H  *  h  >  min  (h  ,h_  )  *  min  (2K_,2a_).  In  this  case 

1 .  St  SS  tt 

* 

Pg  fi  »  0  and  there  is  a  unique  optimal  solution,  x  ,  to 

P_  which  satisfies, 

T 

(i)  x*  €  ?(PgJ.»PCs)  and 

*  * 

(ii)  d(pst,x  )  -  ug  and  d(ptg,x  )  «  . 

* 

That  is,  x  is  unique  and  lies  on  the  path  between  p  and 

st 

p_  a  distance  of  u  units  from  p  or  u  units  from  p 
rts  s  rst  t  rts 

Case  2 :  H  *  h  ■  2j?  . 

-  ss  s 

In  this  case  RC  is  a  connected  sec  (see  Remarks  1  and  2)  which 

is  contained  entirely  in  P  •  More  precisely, 

s 

RC  *  {x  €  p  (max  v  <  d(p  ,x)  <  min  v  }  (11) 

s  t*s  11  "  3  i+s  2i 

where  for  i^s  v^  and  v  are  calculating  by  the  following 
rules,  which  depend  upon  whether  or  not  P^  "  Pg  *  0. 
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Case 

JL tel-  Pi  n  ? 

s  "  01 

then 

vli 

=  max{d(ps,psi) 

•  (ui 

-  D,i>'0] 

V2i 

-  min{d(pg,pgi> 

+  (ni 

-  »,1>.4(P,. 

.qs)}  • 

Case 

-££li  Pt  n  ps 

+  0. 

Let 

?.  Op  *  P(a  . 

>b  .) 

where  a  , 

€  P(p  ,b  .)  and  let 

is  si 

si 

si 

s  si 

vu 

=  max{d (p  ,a  .  ) 
s  si 

•  V 

0} 

V2i 

=  min{d(pg,bgi) 

+  V 

d(ps,qs)}  • 

The  method  is  illustrated  for  Cases  1  and  2  in  the  following  examples 
which  are  taken  from  [  1  ] . 

Examp le  1:  See  Figure  1.  g1  =  6,  g2  *  7,  g3  =*  6,  =»  0, 

°13  =  2’  °23  “  2’  h12  “  81  +  S2  *  13’  ^13  =  2(2  +  6  +  6)/2  ~  14, 

h23  *  2(2  +  7  +  6)/2  ■  15.  Therefore  H  ■  h23  *  15.  This  is  Case  1 

u3  *  (15  -  12)/2  =  1.5,  u2  »  (15  -  14)/2  =  .5  and  x  is  the  unique 

point,  as  shown  in  Figure  1,  .5  units  from  a  and  1.5  units  from  b. 

Exaople  2:  See  Figure  2.  gL  =3,  g2  *  7,  g3  «  5,  D12  *  0, 

°13  *  2’  °23  *  2’  h12  *  *1  +  h  ’  15  ’  hi3  *  2(2  +  8  +  5)/2  -  14, 

h23  *  2^2  +  7  +  5)/2  3  therefore  H  ■  h^  *  2g^  ■  16.  This  is 
Case  2;  u2  »  (16  -  14)/2  -  1,  u3  »  (16  -  10)/2  -  3.  For  i  -  2, 

H  P2  *  P(p2,a)  which  is  Case  2(b),  v^2  *  max{3  -  1,  0}  *  2  and 

v,2  *  min{6  +1,  8}  »  7.  For  i  *  3,  P^  H  P3  *  0.  This  is  case  2(a) 

thus  p^3  *  a  and  v^3  *  max  [6  -  (3  -  2),  9}  *  5, 

v23  -  min( 6  +  (3  -  2),  8}  -  7. 


i 
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RC  »  {x  €  p^:  max  (2,5)  <  d(p^,x)  <  min  (7,7)} 

=  [x  €  p^:  5  <  d(p,x)  <  7],  RC  is  illustrated  in 
in  Figure  2. 

Once  the  distance  matrix  has  been  determined  the  computational 

2 

effort  involved  in  finding  RC  is  0(m  )  operations  to  calculate  H 

plus  whatever  is  necessary  to  calculate  the  v. . ,  a  . ,  and  b  . .  In 

ij  si  si 

the  next  section  we  provide  results  which  justify  an  0(m)  method 
for  locating  an  optimal  solution  to  the  unweighted  version  of  P^. 


3 .  Finding  an  Optimal  Facility  in  an  Unweighted  Tree 

Suppose  that  w(i)  *  c  for  all  i  1  <  i  <  m  or  without  loss  of 
generality  that  c  *  1.  Given  x  €  T  and  v  6  V(x)  recall  the  definition 
of  Ty(x)  from  Section  2  as  the  connected  subtree  of  T  which  contains 
the  segment  (x,v)  plus  all  y  €  I  such  that  v  6  P(y,x).  The  validity 
of  the  following  remark  is  based  on  the  fact  that  T  is  a  tree. 


Remark  3 : 


r  (x)  and  q.  €  T  (x)  if  and  only  if  P.  C  T  (x). 

V  XV  i  —  v 


Definition  1:  L(x)  =*  (x)  , . . . ,  ^  (x) }  is  the  set  of  indices  of 

the  node  pair3  (Pj  (xy  ‘fy,  (x))  where  the  maximum  in  calculating  r(x) 
from  (2)  occurs.  L(x)  represents  the  set  of  bottleneck  node  pairs 

for  each  i,  1  <  i  <  k(x),  is  a 
bottleneck  path  with  respect  to  x. 

Although  we  assume  in  this  section  that  w(i)  *  1  for  all  i, 

1  <  i  <  m,  the  following  lemmas  1  and  2  hold  true  for  any  w(i)  >  0.  We 
state  and  prove  the  lemmas  using  any  positive  weights  since  the  more  general 


with  respect  to  x  and  P^  ^ 


results  will  be  used  later  in  Section  4  in  the  weighted  version  of  PT> 
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Lemma  1:  Given  x  £  T  suppose  that  for  some  l^(x)  €  L(x),  1  <  i  <  k(x), 

and  for  some  v  £  V(x)  Chat  Pi  ,  .  —  T  (x).  Then  RC  c  T  (x). 

•u^(.x;  —  v  —  v 

*  •* 

Proof:  Suppose  there  exists  x  c  RC  -  Tv(x).  Let  r  denote  the  optimal 

k 

round  trip  distance  (i.e.,  r  is  the  optimal  objective  function  value  of  P 
when  w(i)  >  0  for  all  i).  For  notational  convenience  let  s  =  •L^(x).  Since 

•jp 

Pg  C  T^(x)  we  must  have  Pg(x  )  €  (x) ,  and  since  by  assumption  x  j£  T  (x) 

we  have 


(i) 


P(x  ,p  (x  ))  and 


(ii)  p  (x  )  =  p  (x) . 

5  S 


Consequently, 

ie  k  k  ★ 

r(x  )  *  r  >w(s)(d(x  ,  pg)  +  d(x  ,q  )  +  k  ) 

"  2w(s) • (d(x*,pg (x*>)  +•  gg) 

k  k 

3  2w(s )  •  (d(x  ,x)  +  d  (x ,  p  (x  )  +  g  ) 

s  s 

3  2w(s )  •  (d(x  ,x)  +  d(x,pg(x)  +  gg  ) 

* 

*  2w  ( s )  •  d  (x  ,x)  +  r(x) 


k 

by  definition  of  r 

from  (6 ) 

from  (i) 

from  (ii) 

from  (6)  and  s  €  L(x) 
★ 


>  r  (x) 

k 

which  contradicts  the  assumption  that  x  £  RC.  U 


since  x  ^  x 


Corollary:  Given  x  €  T  suppose  that  for  some  vi>v2  ^  V(x),  vj_  *  v  , 

and  some  i,j,  1  <  i  <  j  <  k(x),  that  P.  ,  ,  =  T  (x)  and 
—  —  —  ■v^v,x;~v^ 


PL.(x) 


T  (x).  Then  x  is  the  unique  optimal  solution  to  P_. 
v2 
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Lemma  2 :  Given  x  €  T  and  any  {^(x)  £  L(x)  suppose  Chat  x  €  p^ 
Chen 

(i)  x  £  RC  and 
(li)  «CCP4i(l)  . 


Let 

s  *  -t-i(x).  Consider  any  point  y  +  x. 

We  have 

r  (y)  > 

w(s)rs  (y) 

by  definition  of 

r(y) 

2w(s)  •  (d(y,ps  (y))  +  g^) 

from  (6) 

(12) 

> 

2w(s ) • (d (x, pg (x) )  +  gs ) 

since  p  (x)  »  x 

3 

(13) 

3 

r  (x) 

since  s  €  L(x) . 

This  proves  (i).  From  lines  (12)  and  (13)  it  is  evident  that  if 
d(y,po  (y) )  >  d(x,p,(x))  «  0,  then  r(y)  >  r(x).  Thus  r(y)  ■  r(x) 

only  if  d(y,p  (y))  *0  or  equivalently  only  if  y  €  P  ,  which 

3  3 

proves  (ii).  □ 

Lemma  2  says  that  if  x  €  P^  ^  for  some  V  (x)  €  L(x),  Chen 

x  is  an  optimal  location  and  the  set  of  all  optimal  locations  for  P 
is  concamed  in  the  bottleneck  path  P^ 

Lemma  3 :  r(x)  <  r(y)  +  2d(x,y)  ^  x,y  €  T  . 

Proof:  We  have  from  Che  definition  of  r(y) 

r1(y)  *  d(y,pi)  +  d(y,qi)  +  ^  <  r(y)  V  1  <  t  <  m  (14) 

Also,  making  use  ot  the  triangle  inequality  which  is  valid  for  trees 
we  have  for  all  i, 

d(x,pi)  +  d(x,qt)  +  kt  <  d(x,y)  +  d(y,p1)  +  d(x,y)  +  d(y,qt)  + k  (15) 
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Combining  (14)  and  (15)  we  get  for  all  i, 

dOt.p^  +  dCx.q^  +  ^  <  2d(y,x)  +  r(y)  (16) 

and  finally  since  r(x)  *  max  [d(x,p.)  +  d(x,q.)  +  k.}  we  get 

l<i^n  1  li 

the  desired  result.  □ 

Now  we  would  like  to  define  what  we  shall  call  an  "extreme 
point"  and  an  "interior  point"  of  RC, 

Definition  2:  A  point  x  2  RC  is  an  extreme  point  of  RC 
if  and  only  if  there  exists  a  number  5  >  0,  5  arbitrarily  small, 
and  a  point  y  2  RC  such  that  (i)  d(y,x)  *  5  and 
(ii)  (P(y,x)  -  {*})  0  RC  -  0. 

Remark  4:  From  lemmas  1  and  2  and  their  corallaries  it  is  clear  that 
RC  is  either  a  connected  set  which  is  contained  in  a  bottleneck  path 
or  a  single  point.  In  any  event  RC  can  be  written  as, 

RC  =»  P(a  ,b  )  »  (x|d(a  ,x)  +  d(x,b  )  -  d(a  ,b  )}  (17) 

'k  k 

where  a  and  b  are  the  endpoints  of  RC. 

From  Definition  2  and  Remark  4  it  is  clear  that  if  x  2  E 

Hr  Hr  ^ 

then  x  =  a  ,  x=*b,  or  x  €  V,  which  leads  us  to  the  conclusion 

E  C  (RC  n  V)  U  [a*,b*} . 

Definition  3:  A  point  x  €  I  -  RC  -  E  is  called  an  interior 
point  of  RC. 

Next  we  would  like  to  associate  with  each  point  x  2  T  a 
point  in  RC,  which  we  call  c(x).  We  define  c(x)  as  follows: 


a 
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Definition  4:  Lee  c(x)  be  Che  unique  point  in  RC  which 
is  closest  to  x,  that  is  c(x)  satisfies 

d (x , c (x ) ) - min  [d(x,y)}  , 
y€RC 

Note  that  if  x  €  RC  then  c(x)  -  x. 

Definition  5:  For  all  x  €  T-I  let  v  be  defined  as  any 
■  x 

node  in  V(c(x))  which  satisfies  the  following  two  properties: 

(i)  x  €  Tv  (c(x)) 
x 

(ii)  RC  C  (T  -  Tv  (c  (x ) ) )  U  { c  (x ) }  • 
x 

Later,  in  Lemma  5,  we  show  that  v^  is  well  defined. 

For  an  example  of  the  previous  definitions  see  Figures  1  and  2. 

In  Figure  1  E  -  RC  -  {x  }  and  1-0;  c(x)  ■  x  ^  x  €  T;  for  each 

*  * 
xcT(x),  v  -  a;  for  each  x  €  T,  (x  ) ,  v  -  b. 

SL  X  D  X 

In  Figure  2  RC  is  shown:  E  -  {a  ,a,b  };  for  each  x  5  Tb(a), 

★  * 

c (x)  -a  and  v  *  b;  for  each  x  €  T  (a  ) ,  c (x)  -  a  and 
x  p2 

v  -  p, ;  for  each  x  S  T  (a),  c(x)  -  a  and  v  »  q?;  and  for  each 
★  ★  ^ 

x  c  T  (b  ),  c(x)  -  b  and  v  -  q,  . 
qL  x  1 

In  the  next  lemma  we  provide  an  alternate  characterization  of 

the  extreme  points  of  RC  which  in  turn  leads  to  a  proof  that  for 

all  x  €  T-I,  v^  is  well  defined. 

Lemma  4:  E  -  [x  €  Rc|x  -  c(y)  for  some  y  c  T-l}. 

Proof:  P)  Let  x  -  c(y)  for  some  y  €  T-I.  Now  if  y  €  T  then 
by  definition  x  -  c(y)  -  y  6  E.  Suppose  that  y  t  E.  Then  x  is 
the  closest  point  in  RC  to  y.  Thus  we  have  (i)  d(x,y)  >  0  and 
(ii)  (P(x,y)  -  Cx})  ^  RC  -  0  which  shows  from  Definition  2  that  x  E 
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G>  Let  x  €  E.  From  Definition  2  we  see  that  there  exists 
a  point  y  2  RC  such  that  (P(y,x)  -  {x} )  f!  RC  *  0.  Since  RC  is 
connected  and  T  is  a  tree  this  implies  that  x  is  the  closest  point 
in  RC  to  y.  □ 

Lemma  5:  For  all  x  €  T-I,  v  is  well  defined  (that  is  v  exists). 

x  x 

Proof:  There  are  two  cases  to  consider, 

Case  1:  x  €  E. 

Then  c(x)  ■  x  2  E  and  from  Definition  2  there  exists  y  t  E 
which  is  arbitrarily  close  to  x.  Since  y  is  arbitrarily 
close  there  exists  v  2  V(y)  such  that  v  €  V(x)  and 
y  2  P(v,x).  Clearly  x  6  T^(x)  and  this  proves  (i)  of 
Definition  5.  Next  since  x  €  E  we  have 

(P(y,x)  -  [x])  fl  RC  =  0  and  since  RC  is  connected  we  must 
have  RC  C  (T  -  Tv(x))  U  {x},  which  proves  (ii). 

Case  2:  x  i  E,  and  since  x  €  T-I  =*  x  2  RC 
Since  x  ^  c(x)  and  P(x,c(x))  H  RC  ■  [c(x)}  there  exists 
a  y  2  P(x,c(x))  such  that  y  is  arbitrarily  close  to  c(x). 
We  use  the  same  argument  as  in  Case  1  to  prove  (i)  and  then 
recalling  from  Lemma  4  that  c (x)  €  E  we  again  use  the  same 
argument  as  in  Case  1  to  prove  (ii).  □ 

Remark  5:  Although  Lemma  5  shows  that  v^  exists  for  all  x  €  T-I, 
it  is  not  necessarily  true  that  v^  is  unique.  For  example  consider 
x  *  a  in  Figure  2.  Here  v^  “  ^2  or  va  *  ^  will  satisfy  the 

properties  of  Definition  5.  It  can  be  shown  however  that  if  x  £  RC, 
then  v^  is  unique. 


We  are  now  ready  to  prove  the  following  Important  result  . 
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Lemma  6:  For  all  x  €  E  there  exists  an  i,  1  <  i  <  k(x),  and 
^(x)  -  L(x)  such  that  ^  C  (T-Ty  (x))  'J  {x)  for  all  vx  €  V(x). 

Proof:  Let  s  ■  <«^(x).  Since  x^E  there  exists  y  t  RC  which  is  arbitrarily  close 

to  x  such  that  x  ■  c(y).  Now  suppose  to  the  contrary  that 

for  all  i,  1  <  i  <  k(x),  and  all  v^  €  V(x)  we  have 

Ps  4  (T-Ty  (x))  U  {x}.  Then  we  intend  to  show  that  there 
x 

exists  a  point,  x  .  with  the  property  that  (i)  x  €  RC 
and  (ii)  d(y,xff)  <  d(y,x).  This  of  course  contradicts  the 
fact  that  x  ■  c(y). 

First  without  loss  of  generality  we  assume  that, 

Pg  €  Tv  (x)-Cx]  for  all  i  1  <  t  <  k(x)  (18) 

X 

and  f0r  aii  v  c  V(x). 

X 

Otherwise  Pj  €  (T-Ty  (x))  U  lx}  and  qs  €  (T-Tv  (x))  'J  [x] 

x  vx 

by  Remark  1  we  would  have  Pg  c  (T-Ty  (x))  U  (x}  which  contradicts 

x 

the  supposition  of  the  proof. 

Now  given  any  €  V(x)  we  partition  L(x)  into  two  secs  de¬ 
pending  upon  the  following  two  cases: 

Case  1:  q  €  T  (x)  -  fxl. 

"■  s  v 

X 

Given  (18)  and  Remark  1  we  have 

Pg  S  (x)  ”  [x]  and  thus 
x 

5g  »  d(x,pg(x))  >  0.  (19) 

Let  L. (x)  -  {3  €  L(x) j  5  >  o}. 

*•  s 


-15- 


Case  2:  qs  €  (T-Tv  (x))  U  {x}. 

X 

Given  (18)  this  implies  chat  x  €  P  and  p  (x)  ■  x. 

s  s 

In  this  case  we  lec 

LjCx)  ■  [s  €  L(x)|ia  ■  O}. 

Clearly  L(x)  *  L^(x)  U  L^Cx)  and  L^(x)  fl  L2(x)  -  0. 

Now  we  defiae  several  numbers  which  will  be  used  in  verifying  the  fact 
Chat  xff  has  the  desired  properties  (i)  and  (ii)  stated  previously. 

m  ★  , 

Let  s,  ■  min  i.  (r  -  r  (x))/2j  >  0 
1  IfEL(x)  1 

ie 

i ^  equals  one  half  of  the  minimum  difference  between  r  ,  the  optimal 

round  trip  distance,  and  the  round  trip  distance  from  x  to  all  non- 

bottleneck  node  pairs. 

Lec  e  ■  min  f  5 . }  >0 
2  iSL^x)  L 

equals  the  minimum  distance  between  x  and  any  bottleneck  path 

which  is  completely  contained  in  T  (x)  -  [x] . 

x 

Let  *  min  d(x,p.)i  >0 
i€L2(x)  1 

equals  the  minimum  distance  from  x  to  p^  over  all  i  S  L(x)  such 
that  x  €  p^. 

Now  choose  y  f  RC,  so  chat  y  €  X  (x)  and 

x 

d(y>x)  *  <  min{ ,s ^ ,e ^ »d (v  ,x) } .  Note  that  this  can  be  done  since 

x  €  E . 
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Next  let  e  *  —  and  let  x  be  the  unique  point  In  P(y,x) 
which  is  s  units  from  x. 

Claim:  x  6  RC. 

★ 

We  will  prove  the  claim  by  showing  that  t\  (x„)  <  r  for 
each  i,  1  <  i  <  m. 

Case  1:  i  i  L(x). 

(a)  p t Cacg )  s  Tv  (x)  -  {x}  . 

11  x 

Since  p  (x  )  €  V  and  x  €  P(v  ,x)  it  must  be  true  that 

rl  £  £  X 

Pt(xff)  -  Pt(x)  and  d(p1(xe),x£)  <  d(pi(xe) ,x) ,  which  implies 

★ 

from  (6)  that  r . (x  )  <  r. (x)  <  r  . 

i  £  —  i  — 

(b)  Pi(x{)  €  (T-Tv  (x))  U  [c(x)}  . 

*  x 

Then  x  €  P(xglpt(xg))  which  implies  that  pt(xg)  *  p^x). 

We  have 


r1(xg)  ■  2(d(xg,pt(xg)  +  gt) 

■  2(d(xg,x)  +  d(x,pi(xg))  +  g^) 

»  2(d (xg iX)  +  d(x,pi(x)  +  qt) 

(r  -r . (x) )  r  (x)  * 


from  (6) 
since  x  6  Pte  ,p. (x_>) 

Sir. 

since  Pt(x.)  ■  P^x) 


Case  2 :  i  €  L(x). 

(a)  i  €  L^Cx) 

Since  x„  €  P(x,p.,(x))  we  have  p.  (x)  -  p .  (x  )  and 
£  l  ri  ric 

d(x  . p. (x  )  <  d(x,p  (x) )  and  thus  from  (6)  r, (x  )  <  r . (x)  -  r*. 
cts  i  i  s  i 

(b)  i  €  Lj (x) 

Since  £  <  e,  we  have  x  €  P.  and  thus  r , (x  )  *  r, (x)  *  r  . 
j  el  i  e  i 

Finally  since  c  <  we  have  xg  €  P(y,x).  The  existence  of  xg 

contradicts  the  fact  that  x  S  E.  C 


In  the  next  Lemma  we  establish  a  relationship  between  r(x), 
* 

c (x) >  and  r  for  all  x  €  T. 


Lemma  7 :  For  all  x  €  T  r(x)  -  r  +  2d(x,c(x)). 


Proof:  If  x  €  RC,  Chen  c(x)  ■  x  (See  Definition  4)  and  thus  the 

proposition  is  true.  If  x  i  RC  then  from  Lemma  4  c(x)  ~  E. 

From  Lemma  6  there  exists  an  i,  1  <  i  <  k(c(x)),  ^(cCx)), 

and  v  €  V(c(x))  such  that  P.  ....  —  (T  -  T  (c(x)))  U  fc(x)} 
c(x)  1  c(x) 

and  x  €  Tv  .  (c(x)).  Let  s  -^(cCx)).  Now  since  x  €  Tv  (x)(c(x)) 
c 1  c 

and  p  (c(x))  €  (T  -T  (c(x)))  U  (c(x)}  we  must  have 


(i)  c(x)  €  P(x,p  (c(x)>)  and  (ii)  p  (c(x))  *  p  (x) 

S  S3 


From  this  we  get. 


r  (x)  *  rg(x) 


by  definition  of  r(x) 


2  (d (x,pg (x) ) )  +  gg 


from  (6) 


2  (d  (x , c  (x ) )  +  d(c  (x) ,  pg  (c  (x) ) )  +  gg 


from  (i)  and  (ii) 


2d (x,c (x) )  +  r (c (x) ) 


from  (6)  and  s  €  L(c(x,i) 


■  2d(x,c (x) )  +  r 


since  c (x)  €  RC. 


Combining  this  with  Lemma  3,  setting  y  -  c(x),  we  get 
r (x)  <  r(c(x))  +  2d (x,c (x) )  ■  r*  +  2d(x,c(x))  which  yields  the 


desired  result.  □ 


Leona  7  cells  us  chec  for  any  location  x  -  T,  the  round  crip 
value,  r(x),  of  this  location  is  linearly  proportional  to  the  distance 
of  this  locacion  from  SC.  Thus  all  of  the  points  in  RC  have  the  same 
round  trip  value  and  the  further  we  move  from  RC,  the  larger  the  round 
trip  value  becons.  This  behavior  is  indicative  of  the  fact  that  r(*) 
is  a  comer  function  on  T.  RC  represents  the  only  sec  of  local  and 
global  optimal  solutions.  This  observation  helps  in  the  proof  of  the 
next  lemma. 

Lemma  8 :  For  any  point  x  €  T  and  every  i,  1  <  i  <  k(x), 
c(x)  €  P(x,p^  (x)(x)). 


Proof:  Lee  s  -  4  (c(x)).  If  x€  RC,  then  c(x)  -  x  and  thus  the  proposition  is 

trivially  true.  Now  suppose  that  x  £  RC.  From  Lemma  6  there 

exists  an  s  and  vc(x)  such  chat 

E  *  *  (c(x)l)  U  [c  (x)}  and  x  €  T  (c(x)). 

C(x)  vc(x)  j 

(See  Figure  3).  We  have  at  this  point  (ii)  c(x)  €  P(x,p  (c(x))).  I 

s 


which  satisfies  the  hypothesis  of  the  1 
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★ 

Now  we  use  P  to  show  that  c(x)  €  PCx.p^x))  for  all 
i  €  L(x).  Suppose  to  the  contrary  that  there  exists  a  j  €  L(x) 
such  that  c(x)  2  P(x,p.(x)>.  Then  there  exists  a  node 
q  €  P(x,c(x))  where  P  and  P(x,p^(x))  separate,  (See 
Figure  3).  Then  we  have 


r (x)  -  2(d(x,q)  +  d<q,pj  (x))  +  gj 

since  j  €  L(x) 

(20) 

-  2 (d (x , q )  +  d(q ,c (x) )  +  r* 

from  lemma  7 

(21) 

Rearranging  (20)  and  (21)  we  get 

r  (c  (x) )  =  r  *  2  (d  (q ,  p ..  (x) )  +  -  d(q,c(x))) 

<  2  (d(q  ,  Pj  (x)  )  +  gj  +  d  (q  ,c  (x)  )  ) 

since  d(q,c(x)) 

>  0 

-  r^  (c(x)). 

But  this  contradicts  the  definition  of  r(’)>  □ 

Lemma  8  says  chat  c(x)  is  an  element  of  the  path  from  x  to 
the  path  between  every  bottleneck  node  pair  of  x.  This  property  is 
crucial  in  the  proof  of  the  following  theorem  which  provides  a  simple 
method  for  solving  the  unweighted  version  of  P  . 

Theorem  1: 


Step  1: 

Choose  any 

xx  €T 

Step  2: 

Calculate 

w 

for 

each 

i  1  <  i  <  m 

Choose 

any  t  (x^ ) 

^  L(xx) 

and 

find 

=  x2 

Step  3: 

Calculate 

ri(x2) 

for 

each 

i  1  <  i  <  m 

Choose 

any  (x2 ) 

€  L(*2) 

and 

find 

Vx2>  *  X3 * 

Let  x*  be  the  unique  point  on  P(x2,x3)  such  that 


d(x3’x  )  -  8<-i(x1)  +  d(*2’*3)  '  S<-i(x2) 

- 

Then  calculated  as  in  Steps  1  -  3  is  optimal. 


* 
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Proof:  We  plan  to  show  that  x  €  P^^),  c(x^))  and  since  RC 

*  r 

is  a  connected  set  this  implies  that  x  e  RC. 

First  note  by  construction  that, 


since  x  6 


8<,.(x.)  +  d(x2»x3J  *  8<,.(x.)  +  8t.(x,)  from  Step  3 

L  1 _ L  1  1  2 

2 

+  “(W  -  Hfy) 


*  2 


2(d(x  ,x  )  +  -  r^(Xi> 


(x  ) 


(22) 


Now  from  Lemma  8  we  have  c(x2)  €  P(x2>x3)  and  since  RC 
is  connected  and  T  is  a  tree,  we  have  c(x^)  €  P(x2>x3). 

•ff 

By  construction  x  €  P(x2>x3).  We  claim  that  x  €  P(c  (x2)  ,c  (x^) )  • 
Suppose  to  the  contrary  that  x*  t  P(c (x2 ) ,c (x3) ) . 

Lee  us  txrsc  assume  chat  x*  €  P(c(x2),x2)  -  £ c(x2>} . 

See  Figure  4. 

From  (22)  we  have 


rw 


(x*)  ■  r4  (x  j(x  )  ■  2d(x*,c(x2))  +  r 


from  Lemna  7  (23 ) 

and  since 
★ 

c  (x  )  -  c  (x2 ) 


However  since  x  €  Pfc^)^)  -  £c(x2)3» 

r41(x1)(c(x2))  “  2d(c(x2)’x*)  +  \(Xl 
■  4d  (c  (x2 ) .  x2 )  +  r 


,(x*> 


>  r 


from  (23) 


since  d(c(x2) ,x  )  >  0. 


But  this  contradicts  the  fact  that  c(x2>  €  RC.  The  proof  for  the  case 
*  . 

where  x  c  P(c(x3>,x3)  -  lc(x,)}  is  identical.  □ 


— F 
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Example  3:  We  perform  Che  procedure  of  Che  Cheorem  on  Che 
problems  in  Figures  1  and  2.  First  in  Figure  1: 

Seep  1  Choose  x^  »  p^ 

Step  2  r^p^)  *  12 ;  r2^pl^  *  16;  r3^pl^  =  24; 

Kpj^)  -  {3};  x2  =»  P3  (P 1  >  3  b- 

Step  3  r^(b)  =*  16;  r2(b)  =  18;  r^Cb)  =  12 

L(b)  =  (2}  ;  x3  =■  p2  (b)  =  a 

"it 

x  €  P(a,b),  6+2-7  *  .5  units  from  a. 

2 

In  Figure  2: 

Step  1  Choose  x^  ■  p^ 

Step  2  r^Pj)  3  16»  r2^pl^  =  20;  r3(pl)  -  26 
LCp^  -  C3};  x2  »  P3  (p^)  ■  b 
Step  3  r^b)  =■  20;  r2(b)  *  24;  r3(b)  *  10; 

L(b)  »  {2} ;  x3  »  p2  (b)  *  a. 

x  €  P(a,b),  5+2-7  *  0  units  from  a. 

2 

★ 

The  computational  effort  involved  in  calculating  x  is  0(m) 
which  is  essentially  the  amount  of  work  which  is  required  to  calculate 
r^Cx^)  in  Step  2  and  r^^)  in  Step  3.  In  the  next  section  we 
discuss  a  method  for  solving  the  weighted  version  of  P^,. 

4.  Round  Trip  Location  on  a  Weighted  Tree 

We  assume  that  w(i)  >  0  for  all  i,  1  <  i  <  m.  First  we 
present  a  method  for  locating  a  single  edge  of  T  which  is  guaranteed 
to  contain  an  optimal  solution  to  P^.  Then  we  provide  a  method  for 
finding  a  "local"  solution  to  P  ,  chat  is  a  solution  to  PT  where  a 
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locacion  is  constrained  to  lie  on  a  single  edge  of  the  tree.  Combining 
these  two  approaches  yields  an  algorithm  for  P^. 

4.1  Locating  an  Optimal  Edge 

Using  Lemmas  1  and  2,  which  were  stated  and  proven  for  the 

weighted  version  of  P^,  we  devise  an  algorithm  which  results  in 

exactly  one  of  the  following  states, 

(i)  We  find  a  node  v  6  V  which  satisfies  the  properties  of 

Lemma  2  and  therefore  is  optimal. 

or  (ii)  We  find  a  single  edge  e  =  (v  , v^)  of  T  which  is 

guaranteed  to  contain  an  optimal  solution  to  p  . 

In  case  (ii)  we  show  in  Section  4.2  how  to  find  an  optimal  solution  on 

the  edge  (v  ,v  ). 

st 

First  we  state  the  algorithm  for  (i)  and  (ii),  and  then  we 
justify  its  validity  by  reference  to  Lemmas  1  and  2.  We  assume  that 
T  contains  at  least  one  edge. 

Algorithm  4 . 1 

1.  Set  k  »  0,  To  =  T. 

2.  If  consists  of  only  the  single  edge,  (v^.v^),  Chen  go  to 

5.  Otherwise  choose  v  €  V  such  that  |v(v  )|  >1.  (In 
Algorithm  4.2  co  be  stated  later  we  specify  more  precisely  the 
method  for  choosing  v^. ) 

3.  Calculate  L(vk).  Pick  (^(v^)  £  L(vk),  1  <  i  <  k(vk). 

If  v.  €  p.  .  .  then  go  to  4.  Otherwise  find  v  €  V(v  ) 
i'  k'  K  k 

SUChth4t  *tl(vk)  2  %<V'  L“  IM-IkniVk<V- 

Replace  k  by  k  +  1.  Go  to  2. 

Stop.  v^  £  RC. 

Stop.  RC  ~  (V  ,V  )  *  0. 


5. 


Ac  iteration  k  of  Algorithm  4.1  we  choose  a  node  v^  €  V 

which  is  not  a  "leaf"  of  T.  This  is  so  that  in  Step  3,  T-  (v  )  4  T, 

vk  k 

i.e.  so  that  T_  (v  )  is  a  proper  subtree  of  T.  This  is  essential 
vk  K 

for  the  convergence  of  Algorithm  4.1  to  a  single  edge.  In  Step  3  if 

vk  £  ^  then  by  Lemma  2  v^  £  RC,  and  we  go  to  Step  4.  Otherwise 

i '  k 

from  Lemma  1,  RC  —  T—  (v  )  which  justifies  the  reduction  of  T,  to 
'  Vk  k  k 

Tk+^-  Eventually  we  reach  either  Step  4  or  Step  5. 

The  computational  complexity  of  Algorithm  4.1  is  essentially 

of  0(m)  each  time  Step  3  is  performed,  which  is  the  effort  required 

to  calculate  L(v^)  and  (v^).  Mext  we  present  an  algorithm  of 

k 

complexity  0(n)  for  choosing  v^  which  in  turn  guarantees  that  Step  3 

be  executed  at  most  O(Lgn)  times,  thus  making  the  total  computational 

effort  of  Algorithm  4.1  0 ( (n-Hn)^gn) .  Since  in  most  practical  cases 

m  >  n,  i.e.  the  number  of  node  pairs  tends  to  exceed  the  number  of 

nodes,  the  complexity  would  then  be  0 (m  Lg  m) . 

Definition  6:  For  any  v  £  V  and  v^  £  V(v),  1  <  k  <  ]v(v)| 

let  N(v)  =  max  [ | V  L  x—  (v ) i } .  N(v)  is  the  maximum  number  of 
l<k<',V(v)|  v  k 

if 

nodes  in  any  of  the  subtrees  T—  (v).  A  centroid .  v  ,  of  T  is  a 


node  with  the  property  that. 


N(v  )  “  ®in{N(v)} 

v€v 


(24) 


*  r  n*| 

Clearly  N(v  )  <  '  +  1,  thus  if  at  each  iteration  of  Step 

L“'  .  flTk  P.  Vi] 

!  V|  <  [-V- 


we  choose  a  centroid  of  T^,  then  we  have 


2 

+1. 


This  implies  that  T.  will  be  reduced  to  an  edge  of  T  in  at  most 


will 


O(Lgn)  operations,  so  that  Step  3  will  be  executed  at  most  O(Lgn)  times. 


-25- 


r(xc*(e))  * 

min 

0<t<d 

(r  (xt(e))} 

(25) 

How  first  we  note 

that 

E^i  (xt  (®  )  ) 

can  be  redefined  from  (6)  as 

r.(xc(e))  »  min{min{ t  +  d  (va  ,  p^  ,d-t-l-d(vt ,  p^ } 

-Hnin{t  +  d(v  ,q  )  ,d-t+d (v  ,q  )}} 

Si  Cl 

+  8t  (26) 

The  following  lemma  uses  (26)  Co  characterize  r^(-)- 

Lemma  9:  For  each  i,  1  <  i  <  m,  r^x^e))  is  linear  in  c, 

0  <  t  <  d,  with  slope  belonging  to  the  set  [2,  0,  -2}. 

Proof : 

Case  1:  xt(e)  € 

Then  either  v  €  p(x  (e),p  )  or  v  €  P(x„  (e)  ,q  . ) . 

5  LI  S  t  1 

Without  loss  of  generality  suppose  that  the  former  case  is  true. 
Then  S  P(x(.(e),qi)  and  from  (26) 

Ci (xt (e) )  *  C+d(vs>pt)  +  d-t+d(vt>qi)  +  gi 
*  d(p.,q.)  +  gt  . 

Thus  r^(*)  is  linear  in  t  with  zero  slope. 

Case  2:  xt(e)  j£  P^ 

Then  since  T  is  a  tree  either  (i)  v  €  P(x  (e),p.(x  (e))) 

s  c  i  c 

or  (ii)  c  P(x(;(e),pi(xc(e))).  If  (i)  holds  true  Chen 

ri  (xt  (e) )  "  c  +  d(v3>Pi)  +  c  +  d(v3'<li)  +  *i 

■  2t  +  d (v  ,p  )  +  d (v  ,q  )  +  q. 

SI  SI  1 

and  thus  r^-)  is  linear  with  slope  2.  In  Case  (ii) 

ri^xt^e^  *  +  d^vt’^i^  +  d(vt>p^)  +  8^  which  is  linear 

with  slope  -2.  O 
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Mow  we  define 

R . (t )  *  max  (w(j)-r . (x  (e))}  .  (27) 

1  l<i<i  J 

Moce  that  R  (t)  ®  r(x_(e))  and  thus  from  (25)  our  objective 
01  C 

is  to  find  the  minimum  of  R  (t)  on  the  interval  (0,d). 

tu 

Corollary:  For  each  i,  1  <  i  <  m,  R^t)  is  a  piecewise  linear 
convex  function  with  slopes  belonging  to  the  set 
[+2w ( 1 ) , . . . , +2w (i ) , 0}  . 

Proof:  From  Lemma  9  each  r  (•)  is  linear  and  therefore  convex. 

Thus  w(j)rj(-)  is  also  linear,  convex,  and  has  a  slope  of 

+2w(j)  or  0.  Finally  it  is  well  known  that  the  maximum  of 
a  set  of  convex  functions  is  convex.  □ 


Remark  6:  Since  R  (t)  is  piecewise  linear  and  convex  we  could  find 

. 1  111  m 

its  minimum  by  solving  a  linear  programming  problem  which  is  equivalent 
to  a  two  person  zero  sum  bi-matrix  game.  For  the  details  see  [4] . 

Instead  of  solving  a  linear  program  we  will  present  an  iterative  method 
which  generates  Rm(t)  and  its  minimum  values  in  0(m  -Lg  m)  operations. 

Next  we  outline  an  interative  method  which  can  be  used  to  generate 
R^+^  from  R^(t).  Since  R^(t)  is  a  piecewise  linear  convex  function 
let  T.  »  i  0  »  t*1,  tj,,...,c*'  -  d]  be  the  set  of  points  in  the  interval 

1  0  1  Q . 

1 

(0,d)  where  the  function  R^(t)  changes  slope.  Then  for  each  k, 

1  <  k  <  ni  and  each  t  €  (t^^.t^)  ve  have, 


Ri(t)  -  VW  +  W  -  R(tk-1}  i  . 

i  i  '  ck-l; 

Ck  *  Vl 
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The  function  R^(t)  is  completely  described  by  the  list  and  the 

values  of  R^(0  at  the  points  of  T^.  Next,  given  the  list  and 
the  values  of  each  member  of  we  outline  how  R  ( •)  can  be 

generated  from  R^O)  or  effectively  how  to  generate  T  ^  and  its 
values  from  T^. 

From  (28)  Ri+^(t)  *  max(Ri(t),w(i+l)r  ^(x  (e))}.  Now  since 

w(i+l)  >...>w(l)  it  must  be  true  that  R^O)  and  w(i+l) -r  ^ (• ) 
either  coincide  with  one  another  or  they  intersect  in  at  most  two 
points.  Figures  5,  6,  and  7  characterize  the  possibilities  which  can 
occur.  In  (i)-(iv)  below  we  describe  how  can  be  generated  from 

T.,  by  taking  into  consideration  the  various  possibilities  which  can 


Let  t(t)  €  [o,d]  be  a  point  where  w(i-fl)  reaches  its 

minimum  (maximum)  value.  Also  let  7  £  T,  be  a  point  where  R^(-) 
reaches  its  minimum  value,  let  Ti/Ci,l:2^  be  the  3et  ail  points  in 
between  t^  and  t We  break  the  different  possibilities  into 
four  cases : 

(i)  If  w(i+l)r^+^(x- (r))  >  R^(t)  (Figures  5 (a) , 6 (a) , 7  (a) )  then 

Vl(,)  =  ri+l(°  and  Ti+I  =  (0>d}' 

(ii)  If  w(i+l)r  (x*(e))  <  R.(t)  and  w(i+l)r  ,  (x-(e))  <  R.  (t) 

x+i  t  —i  l+i  t  — ■  r 

(Figures  5  (c ) , 6 (c ) , 7  (c) )  then 

W°  *  V)  and  Ti+i  -  V 

(iii)  If  w(i+l)ri+1(x^(e))  >  R^t)  and  ri+1(x-(e))  <Ri+i^) 

(Figures  6(b)  and  7(b))  then  we  find  the  points  t^_1  €  x 

i  *  i  i 

and  t^  €  and  the  unique  point  t  ^  ^k-l^k^  suc^  that 


‘  i+1 '  t 


(x  *(e) )  ■  R. (t  ).  Then 


.  •  -  ;*  ■’ 1 

i  *aXts 
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T 


i+1 


CO, t*}  U  T.(t\d)  if  slope  ■  -2w(i+l) 
1 

Ti(0,tk-1)  U  £c*>d}  if  sl°Pe  ”  2w(t+1> 


(iv)  ri+1(x»(e))  *  ri+^  (x^  (e) ) ;  Ri(t)  <  w(i+l)ri+1(x£(e)) ; 


and  w(i+l)r1+1(xQ(e))  <  Ri(0)  or  w(i+l)r1+1 (xd(e) )  <  R^d)  (Figure  5(b)). 


Then  che 
i  * 


etcher  (a)  intersect  at  one  point,  t  , 


t^j.  <  t  <  t^,  or  (b)  they  intersect  at  two  points  t^  and  t2 


where 

★ 

C1  < 

★ 

t2  and 

Ck-1  -  C 

k 

i* 

i 

Ck 

and 

i 

Ci-1 

< 

* 

C2  * 

t1.  Then 
J 

(  T. 

i 

(Q.if) 

U  {c*,d} 

if 

(a) 

and 

ri+l 

(xo 

(•)) 

<  Ri(0) 

to 

,t*}  U  T, 

if 

(a) 

and 

ri+l 

(xd 

(e)) 

<  R1(d) 

T  -< 

i+1 

T. 

i 

<°«4-i) 

U  tc*,c*i 

U  I 

'i<ci 

>  d) 

if 

(b) 

• 

\ 

In  cases  (i)  -  (iv)  a  simple  updating  procedure  has  been  out¬ 
lined  to  generate  Ri+^0)  from  R^t1)*  The  process  is  iterative,  start¬ 
ing  with  R, ( • )  ■  r  (•)  and  terminating  with  R  (•)•  The  interval  ob- 
11  m 

tained  from  che  break  points  in  T  which  minimizes  R  (•)  provides 

m  m 

a  sec  of  optimal  solutions  to  P^.  Notice  that  in  cases  (i)  -  (iv) 

at  most  two  ooincs  are  added  to  T.  at  any  iteration  and  thus  T 

1  cn 

contains  at  most  2m  points.  Also  the  lists  T^  and  their  values 
can  be  stored  in  a  data  structure  which  is  an  ordered  tree  of  depth 
at  most  0(-tg  m).  Thus  we  require  at  most  0(4>g  m)  operations  to  retrieve 
any  information  which  is  necessary  in  (i)  -  (iv).  This  makes  compu¬ 
tational  effort  involved  in  the  iterative  application  of  (i)  -  (iv)  at 
most  0(m  ig  m) .  Combining  this  with  Algorithm  4.1  yields  an  algorithm 
for  PT  of  0((n+mKgn  +  rotgm)  of  if  m  >  n  of  0(m  tg  m). 


Example  4:  We  perform  che  procedure  outlined  in  Section  4 

on  the  problem  in  Figure  8  which  is  identical  to  Che  problem  in  Figure 
2  except  chat  positive  node  weights  have  been  added.  We  describe  the 
calculations  in  Algorithm  4.1. 


Step 

1. 

tQ  -  T 

,  k  -  0 

Step 

2. 

Choose 

node 

a  in  Algorithm  4.2  since  N(a)  -  4. 

Step 

3. 

r^(a) 

-2(6 

+  2 

+  8)  -  24 

r2(a) 

-  3(3 

+  4 

+  7)  -  52 

r3(a) 

-4(3 

+  6 

+  5)  -  56 

so 

L(a)  - 

{3} 

and 

a  j£  P  .  Tj_  -  Tb(a). 

Step 

2. 

Choose 

node 

b 

since  N(b)  -  2  in  Tb(a). 

Step 

3. 

r^(b) 

-2(8 

+  4 

+  8)  -  40 

r2(b) 

-  3(5 

+  6 

+  7)  -  54 

r3(b) 

-4(1+4 

+5)  *  40 

L(b)  - 

•  (2} 

and 

b  jE  P, .  T,  *  T  (a)  H  T  (b)  -  (a,b) 
2  2b  a 

Step 

2. 

t2  -  (v3,vt 

)  - 

(a,b) . 

Now  we  have  limited  RC  to  the  edge  (a,b).  Next  we  apply 
(i)  -  (iv)  to  calculate  R3(-).  R^t)  -  r^x^e))  -  2t  +  32,  from 
(27)  and  *  C 0 , 2 } .  See  Figure  9.  Next  r2(xt(e))  «  3t  +  52  and 
this  is  case  (i)  and  thus  -  Co , 2 } .  Finally  r3(xc(e))  -  -4t  +  56. 

this  is  case  (iii),  t*  -  J  and  *3  -  {0,y,2j.  *3<*)  is  shown  in 

ic  4 

Figure  9.  *  is  located  on  edge  (a,b),  j  units  from  a. 


5 .  Concluding  Remarks 

In  this  paper  ve  have  shown  how  the  special  structure  of  a  tree 
network  can  be  exploited  to  provide  simple  procedures  for  solving  the 
unweighted  and  weighted  round  trip  location  problems.  This  includes. 


as  a  special  case,  Che  one  center  problem  on  a  tree.  In  fact,  the 
method  proposed  in  Section  3  for  solving  the  unweighted  version  of 
reduces  to  the  method  proposed  by  Handler  for  the  unweighted  one  center 
problem  when  pi  *  qi  and  ^  -  0  for  all  i,  1  <  i  <  m.  Although 
the  results  in  this  paper  pertain  to  the  round  trip  location  problem, 
we  feel  certain  that  many  of  the  properties  developed  here  will  hold  true 
for  several  related  problems.  In  particular  consider  a  vehicle  routing 
problem  where  a  vehicle  travels  from  a  facility  location  to  a  set,  S^, 
of  neighborhoods  (or  nodes),  visiting  each  neighborhood  in  S.  exactly 
once  and  then  returning  to  the  facility.  In  this  paper  but 

in  general  S^  might  contain  more  than  two  neighborhoods.  Then  one  might 
consider  locating  a  facility  so  as  to  minimize  the  total  cost  of  visiting 
all  of  the  neighborhood  sets  or  to  minimize  the  maximum  cost  of  visiting  any 
neighborhood  set.  To  our  knowledge  very  little  work  has  been  done  on  combined 
vehicle  routing- location  problems. 

Finally  we  would  like  to  point  out  what  we  consider  to  be  the 
crucial  property  of  the  unweighted  version  of  which  justifies  the 

0(m)  algorithm  described  in  Section  4.3.  Lemma  8  of  Section  2  says 
that  for  each  x  €  T  the  set  of  optimal  solutions  to  PT  intersects 
the  path  from  x  to  each  of  its  bottleneck  node  pairs.  This  property 
no  longer  holds  true  for  the  weighted  version  of  P^.  Consider  the 
simple  three  node  example  in  Figure  10.  Here  p^  *  and  k^  -  0 
for  all  i.  Consider  node  1,  r^d)  *  2  and  r^d)  *  4;  thus 
r(l)  *  100*2  ■  200  and  L (1 )  ■  {2}.  However  the  .vnique  optimal 
solution,  x*.  is  located  between  nodes  2  and  3,  .2  units  from 


node  2. 
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